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ABSTRACT 

We investigate the origin of the galaxy color-concentration bimodality at the bright-end of the 
luminosity function {Mb — 5 log /i 70 < — 18 mag) with regard to the bulge-disc nature of galaxies. Via 
(2D) surface brightness profile modeling with GIM2D, we subdivide the local galaxy population in the 
Millennium Galaxy Catalogue into one-component and two-component systems. We reveal that one- 
component (elliptical and disc-only) systems define the two peaks of the galaxy color-concentration 
distribution (with total stellar mass densities of 0.7 ± 0.1 and 1.3 ± 0.1 x 10^ hjo Mq Mpc~^ 
respectively), while two-component systems contribute to both a bridging population and the red, 
concentrated peak (with total stellar mass densities of 1.1 ± 0.1 and 1.8 ± 0.2 x 10* /170 Mq Mpc~^ 
respectively). Moreover, luminous, 'bulge-less, red discs' and 'disc-less, blue bulges' (blue eUipticals) 
are exceptionally rare (with wo^itme-densities of 1.7 ± 0.3 and 1.1 ± 0.1 x 10~^ /lyg Mpc~^ respec- 
tively). Finally, within the two-component population we confirm a previously-reported correlation 
between bulge and disc color (with a mean offset of only {{u — r)buigc — (w — ?')disc) — 0.22±0.02 mag). 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: fundamental parameters 



1. INTRODUCTION 

The rigorous, statistical confirmation of bimodality 
within the luminous {Mb — 5 log /i7o < — 18 mag) , nearby 
galaxy population has proven a key result to emerge 
from modern wide-field imaging surveys — implying the 
existence of two principal evolutionary pathways for 
bright systems. This bimodality is clearly evident 
in the distribution of local galaxies in various global 

i i.e., total) obs ervable pairs, including color-magnitude 
Strateva et al. 2001; BaMry et a l- .200 4 .), conc entration 

et al. 20Q1 
(Driver et all 



(log n)-magnitude5 (jPriver et al.ll2006i : Bal 
and global color-concentration (log 



|2006[) . Importantly, many galaxies also display distinct 
identifiable structural components — frequently (in mas- 
sive systems) a central bulge and a. surrounding disc 
(|de VaucouleursI [l959l : lFreemanlll970D . as well as bars, 
nuclei, and spiral arms. 

Hierarchical c luster ing (iWhite fc Reed 
White fc Fren kI [l99H iKauff mann et al ' 



1978 



1999 



Somcrville & Primackl Ii999[ ). the predominant the 
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total galaxy light profile. In the context of our two-component 
model, where only the central excess above an outer exponential 
(i.e., the bulge) is fit with a Sersic function, the corresponding 
index will be denoted rit,. 



ory of galaxy formation, postulates the existence of two 
distinct structural formation mechanisms — the cooling 
of gas inside rotating da rk matter halos to form discs 
(jFall fc Efstathioul I198O0 and the merging of similar- 
sized dis cs to form spheroids (i.e., ellipticals and classical 
bulges) (T oomre fc Toomrel[l97a iBaugh. Cole, fc FrenkI 
1996). Hence, understanding the role of structure in 
shaping the observed global bimodality may prove 
crucial to understanding a range of galaxy formation 
processes. 

Using th e Millennium Galaxy Catalogue (MGC, 
iLiske et al. 2003) of 10,095 galaxies. Driver et al. (20061) 
revealed a clear relationship between visual morphologi- 
cal type and position in the global color-concentration 
plane. Specifically, E/SOs dominate the red, high- 
n (centrally-concentrated) peak, Sd/Irrs dominate the 
blue, low-n (diffuse) peak, and Sa-Sc galaxy types span 
both. The authors proposed that this relationship de- 
rives from the two-component nature of galaxies, and 
that objects falling between the two peaks of the global 
bimodality result from the mixing of separate red, com- 
pact, bulg e and blue, diff u se, di sc components in varying 
degrees. Drorv fc Fisher! (j2007t ) challenged this hypoth- 
esis based upon a sample of 39 SO-Sc galaxies with high- 
resolution, HST imaging. Specifically, they recovered a 
correlation between galaxy color and (visually-classified) 
bulge type, such that classical bulge systems are globally 
red (i.e., also have red discs) and pseudobulge systems 
are globally blue, regardless of the bulge-to-total flux ra- 
tio. 

Here we further explore the origin of the color- 
concentration bimodality using the results of Allen et 
al.'s (2006) bulge-disc decomposition of 10,095 bright 
galaxies in the MGC. In Section [2] we review the 
MGC structural dataset and describe the construction 
of a robust sample, supplemented with rest-frame 
component colors. In Section [3] we contrast the global 
color-concentration distributions of one-component 
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and two-component galaxies; identify key structural 
types and construct their luminosity functions; and 
examine the relationship between bulge and disc colors. 
Finally, in Section [4] we summarize our conclusions, 
discuss our results in light of contemporary galaxy 
formation theory, and reflect upon the challenges facing 
low resolution, automated structural decomposition 
studies. Unless otherwise stated, all magnitudes are 
given in the AB system and a cosmological model with 
= 0.3, nA = 0.7 and Hq = 70 km s^^ Mpc'^ is used 
throughout. 

2. DATA: THE MILLENNIUM GALAXY CATALOGUE 

The Millennium Galaxy Catalogue (MGC) is a 
deep (/uiim,B = 26 mag arcsec"^), wide-field (30.88 
deg^), B-band imaging survey obtained with the Wide 
Field Camera on the 2.5-m Isaac Newton Telescope 
(|Liske et al.l 120031). It i s fully contained within both 
the 2dFGRS (iColless et al. 200l) and SDSS (DR5, 
lAdelman-McCarthv et al.. ,2007 ) , which provide photo- 
metric information in additional filters and ~4,700 high 
quality, spectroscopic redshifts. A follow up campaign 
of sp ectroscopic observations (described in Driver et 
120051) has ensured over 96% redshift completeness for 
a sub-sample (MGC-BRIGHT) of 10,095 galaxies with 
B < 20 mag. 

2.1. Bulge-Disc Decomposition 

lAUen et all ()20060 modeled the 2D light distributions 
of the entire MGC-BRIGHT sample, providing a publicly 
available structural catalogue^. The 2D light distribution 
of each galaxy was fit with both a global Sersic (one- 
component) model and a Sersic bulge plus exponential 
disc (two -component) model^ using t he GIM2D analysis 
package (|Simardlll99l . In addition. [Allen et all (|2006f ) 
designed a logical filter to identify and replace illogical 
fits (e.g. inverted profiles where the inner bulge was fit 
by the exponential component and the outer disc by the 
Sersic component). Based on the best-fit model param- 
eters, the logical filter also offers an objective and quan- 
titative subdivision of the local galaxy population into 
three morphological classes. Namely, 'bulge-plus-disc', 
for which the two-component model is preferred; and 
'elliptical' and 'disc-only', for which the one-component 
model is preferred. The distinction between the latter 
two classes is simply made according to the value of the 
best-fit global Sersic index (i.e., greater than or less than 
1.5 respectively), recalling the absence of dwarf galaxies 
in MGC-BRIGHT. 

2.1.1. Treatment of Elliptical Galaxies 

A known problem exists for automated photometric de- 
composition studies concerning the treatment of elliptical 

® http:/ /www.eso.org/r^jliskc/mgc/ 
Specifically, the Sersic function describes a major axis intensity 
profile of the form: IbiR) = h exp [(iJ/iJe)^/" - l] } where 

the effective radius, _Re, encloses half the component luminosity, 
and /e is the intensity at the effective radius. The exponential 
disc model may be considered a specific case of the Sersic profile 
for n = 1 with the following major axis intensity profile: Id{R) = 
Iq exp{—R/h) where Ig is the central intensity and h = 1.678-Re 
is the disc scale-length. For a com prehensive review of the Sersic 
profile see lGraham fc Driveri II2005I') . 



galaxies. Namely, that two-component fits to eyeball- 
morphological ellipticals will readily converge upon a so- 
lution containing an outer disc, although no such struc- 
tural component is actually present (see iTruiillo et al.1 
12004 iGutierrez et al.l [2003 : 1 Allen et all 12006V Rather, 
these 'false discs' serve only to correct slight deviations 
in the measured light profile from a pure global Sersic 
model. The origin of these deviations may be real (e.g. 
a recent merger event, twisting of the inner isophotes, or 
the presence of an outer halo) or artificial (e.g. incorrect 
e stimation of t he sk y background). 

lAUen et al.l (|2006l ) attempt to remove 'false disc' sys- 
tems from their two-component population via the logi- 
cal filter, which assigns all galaxies with model fits corre- 
sponding to bulge-to-total flux ratios (B/T) greater than 
0.8 to the one-component population. The authors claim 
to demonstrate the success of this approach by contrast- 
ing the inclination distribution of their real disc sample, 
which is roughly consistent with a random distribution 
of projections on the sky, against that of their false disc 
sample, which is biased towards a more face-on distribu- 
tion. 

How ever, the B/T > 0.8 cut employed bv I Allen et al.l 
(|200a is significantly less restrictive than the B/T > 0.6 
cut more c ommonly employed in contemporary studies 
(e.g. iTrujil lo et al. 2004; Gutierrez et al. 2004). Hence, 
we conducted visual inspection of the real (MGC) and 
model (GIM2D), _B-band images and surface brightness 
profiles of 100 bright {B < 19 mag) galaxies from Allen 
et al.'s (2006) bulgc-plus-disc sample in order to gauge 
the reliability of this dataset. A total of 23 probable 
'false discs' were thereby identified amongst 72 real discs 
and 5 systems of indeterminate type. Three examples 
of these 'false disc' candidates are displayed in Fig. [6] 
in the Appendix. These galaxies display elliptical visual 
morphologies, although slight flux asymmetries are re- 
vealed upon inspection of the image residuals from their 
one-component fits. Such asymmetries may derive from 
a diverse range of factors, including recent mergers and 
the presence of undigested satellites. 

Reviewing the inclination distribution of Allen et al.'s 
(2006) bulge-plus-disc sample (their Fig. 9), we note that 
the logical filter has reduced, but not entirely removed, 
the bias towards face-on 'disc' fits. Interestingly, this 
bias can be further reduced by applying the more con- 
ventional B/T > 0.6 cut (as shown in our Fig. [T]). Of 
the 100 visually-classified systems discussed above, 15 of 
the 23 'false disc' candidates and only 2 of the 72 real 
discs are removed by this cut — revealing this to be a ro- 
bust, yet conservative, step towards reducing 'false disc' 
contamination. All galaxies from Allen et al.'s (2006) 
bulge-plus-disc sample with B/T > 0.6 (581 of 4427 ob- 
jects) were therefore reassigned to the one component 
population for the remainder of this investigation. 

2.2. Sample Selection 

Motivated by Driver et al.'s (2005) study of the MGC 
luminosity-size plane, we select only those galaxies with 
(i) 0.013 < z < 0.18 to ensure minimal impact from pe- 
culiar velocities, galaxy evolution, and the K-correction, 
and (ii) (^)g g < 25.25 mag arcsec"^ to exclude the 
lowest signal-to-noise systems for which the photome- 
try is less reliable. We also restrict our analysis to near 
face-on discs (i < 60 deg) to minimise potential bi- 
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Fig. 1. — The inclination distribution of galaxy 'discs' in Allen 
et al.'s (2006) bulge-plus-disc class for systems with B/T < 0.6 
{top panel) and 0.6 < B/T < 0.8 {bottom panel). The former 
is consistent with a random distribution of projections on the 
sky, while the latter is biased towards a more 'face-on' distribu- 
tion, indicative of 'false disc' fits. Postage-stamp, MGC B-band 
images of example galaxies from each B/T interval — spanning a 
range of inclinations — are also presented for comparison. The lower 
B/T fiux ratio systems display lenticular and spiral morphologies, 
whereas the higher B/T flux ratio systems generally display ellip- 
tical morphologies. 

ases due to internal, inclination dependent dust extinc- 
tion (jTuffs et al.ll2004t iDriver et al.ll2007bl ). A limit of 
Mb — 51og/i7o < — 18 mag is then used to isolate the 
well-sampled bright-end of the luminosity function. 

By duplicating their bulge-disc decomposition analysis 
on a sub-sample of 682 galaxies with repeat observations, 
lAllen et all ()2006[ ) establish limits on the robustness of 
their derived global structural parameters against varia- 
tions in the local observing conditions. Namely, for the 
accurate recovery of global properties R,, > 0.5T is re- 
quired (where Re is the global half light radius and F 
the characteristic seeing FWHM); and for the accurate 
recovery of bulge properties i?e,& > 0.5r is also required 
(where i?e,b is the bulge half light radius). Hence, we 
apply the former limit to our entire sample used to con- 
struct global color-concentration distributions, and we 
employ the latter limit to identify bulge-plus-disc sys- 
tems with 'confident' bulge Sersic indices (and hence B/T 
flux ratios) during the computation of component colors 
(see Section [23)1 . 

In our final sample there are 917 elliptical, 1544 disc- 
only, and 1526 bulge-plus-disc systems (812 with 'con- 
fident' bulge Sersic indices). Vi sibility- volume weights 
(jPhillips. Davies. fc Disnevlll99Cl| ) were computed to ac- 
count for the relevant observational biases against each 
structural type. These weights are used throughout this 
paper to recover volume-density distributions from the 
raw data. 

2.3. Global and Component Colors 



All galaxies in MGC-BRIGHT have well- 
matched counterparts in the SDSS-DR5 
(lAdil man-McCarthv et all I2007D . which provides 
PSF, Petrosian, and model magnitud es in the Sloan 
ugriz filter system (|Gunn et al.l Il998f ). We compute 
(apparent) global (u — r)g colors for our entire sample 
using the SDSS model magnitudes because these offer 
the highest signal-to-noise me asurements, p articularly 
for the problematic u-band (seel Baldrv et al.l [2004'). 

We also derive (apparent) component {u — r) colors for 
our bulge-plus-disc systems as follows. Core colors are 
computed from the SDSS PSF magnitudes via the re- 
lation [u - r)core = wpsF " ''PSF + 5 log^Q (Fu/Fj.) . The 
third term in this equation, 5 log2o(F„/Fr), is required in 
order to account for an intrinsic difference in area (and 
hence, enclosed core flux) between the SDSS u and r- 
band PSFs. We then equate this core color with the color 
of the bulge, which is a reasonable assumption provided 
the bulge is a similar size to the SDSS u and r-band PSFs 
(F„ « 1.6 and F^ « 1.4 arcsec respectively). For bulges 
much smaller than the SDSS PSFs the core color may be 
biased by the inclusion of significant disc flux. For bulges 
much larger than the SDSS PSFs use of the core color 
may give undue weight to small-scale, radial variations 
in the bulge color resulting from compact central nu- 
clei, s tar- forming rings, nucle ar bars, and/or nuclear dust 
lanes (iJoeee Kennev. fc Smit h 1999; Carollo et al.ll2001l : 
lFisheijl2006i r The radius inside which the model bulge 
intensity exceeds that of the model disc (i.e., the 'anchor 
point', Tanc) providcs a useful measure of bulge size in 
this context. Of our 812 component bulges with 'confi- 
dent' Sersic indices, 145 are a similar size to the SDSS u 
and r-band PSFs (i.e, 0.75 < ranc/ (Fu,rr) < 1.5), while 
217 are smaller and 450 are larger. Hence, a distinction 
is made in our later analysis between galaxies with 'con- 
fident' and 'less confident' bulge (and disc) colors based 
on this size criterion. 

We then suppose that the (unknown) g-band B /T flux 
ratios of our two-component systems may be approxi- 
mated by their (known) i3-band values, given the over- 
lapping wavelength ranges covered by these two filters. 
This enables the computation of g-band bulge magni- 
tudes from the SDSS g-band model magnitudes, and thus 
u and r-band bulge magnitudes using our bulge colors. 
The corresponding disc magnitudes in these Sloan filters 
are then obtained by subtracting the bulge flux from to- 
tal (model) flux, and disc colors derived from their dif- 
ference. The relevant formula is 




(1) 



-0.4it„ 



S2^Q-0.4(g„odrf-(g-ti)i,) 



X0-0.4r,„odcl _ ^10-0.4(g,„odol-(3-r)b) 



K- and e- corrections are required for the computation 
of absolute magnitudes and rest-frame colors. Individual 
K-correcti ons for galax i es in M GC-BRIGHT were com- 
puted by iDriver et al.l (|2005D by matching broadband 
SEDs provided by the (global) MGC and S PSS mag- 
nitude s to the 27 spectral templates given in iPoggiantil 
(|1997f ). Here we repeat this process in order to generate 
accurate and consistent K-corrections for the bulges and 
discs of our two-component systems using t he separate 
colors derived above. For the e-correction iDriver et al.l 
(|2005f) adopt a pure luminosity evolution model of the 
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form 



(2) 



where (3 = 0.75 in the B-band, regardless of galaxy type. 
For consistency with earlier work we continue to utilize 
the MGC S-band absolute magnitudes derived in this 
manner, but extend this approach to better estimate rest- 
frame {u — r) colors. Specificially, we adjust the (3 values 
for each galax y and cornponen t according to its best-fit 
spectral type. Poggianti| (|1997f ) employ three stellar pop- 
ulation models, roughly corresponding to E, Sa and Sc 
spectral types, as the basis of their 27 templates. The 
low redshift evolution in these models may be appoxi- 
mated by j3 values of 1.1, 1.4, and 2.2 in the u-band and 
0.82, 0.92, and 1.1 in the r-band respectively. 

2.4. Dust 

Accounting for the impact of both internal and exter- 
nal dust extinction is an essential step towards the recov- 
ery of colors and luminosities representativ e of galaxy 
stellar populations. All MGC magnitudes (' Liske et al.l 
[2003.) are corrected for external dust extinction (i.e., that 
ar ising due to dust within the Milky Wa y) using the maps 
of lSchlegel. Finkbeiner. fc Davis! pM) . 

Recentlv. [Driver et al.l (|2007t ) recovered empirical 
inclination-attentuation relations from the MGC struc- 
tural catalogue, and used these to calibrate opacities 
in Tuffs et al.'s (2004) multi-component dust models. 
The mean dust c orrections thereby derived (see Table 
1 of [Driver et al.l [200^ are employed in this study to 
account for the effects of internal dust extinction on 
various disc galaxy observables. For our inclination- 
limited sample {i < 60 deg) we obtain median correc- 
tions of A{u — r)d = —0.28 mag and AMb = —0.27 mag 
for component discs, and A{u — r)b = —0.20 mag and 
AMb — —0.98 for component bulges. These corrections 
are applied to all our disc-only and bulge-plus-disc sys- 
tems, and their impact is indicated via reddening vectors 
marked on the relevant plots in Section [3] 

However, observations indicating a possible varia- 
tion in dust content with morphological type raise 
questions regarding the validity of these corrections 
for the entire bulge-plus-disc galaxy population. In 
particular, a strong decrease is observed in the visually- 
estimated dust content of galaxies from late to early 
Hubble types ([van den Bergh[ [2007f) . although far-IR 
emissio n estimates of dust content su ggest a weaker 
trend ([Vlahakis. Dunne, k Ealei [2005f l. As Tuffs et 
al.'s (2004) dust model was devised with reference 
to nearby spiral galaxies only, and Driver et al.'s 
(2007b) inclination attenuation relation is dominated 
by late-type discs, it may not be appropriate to apply 
the resulting corrections to early-type discs. Although, 
some recent observations suggest a surprising level of 
similarity between the distributions of dust in these 
morphological type s, including the presence of 'dust 
spiral arms' in SOs ([Pahre et al.|[2004 ). In Section [3| we 
discuss this potential source of bias in our results, and 
compute alternative luminosity functions for early-type 
bulge-plus-disc systems. 



2.5. Measurement Uncertainties 



Allen et al.'s (2006) catalogue includes 68% confidence 
limits on the measurement of each structural parame- 
ter, estimated by GIM2D based on the shape of the 
minimum in which each fit converges. The median er- 
ror on the global Sersic indices of galaxies in our sam- 
ple is Alogn — ±0.12. Measurement uncertainties on 
our absolute {u — r) colors may be estimated as follows. 
Firstly, we note that for all galaxies in our sample the rel- 
evant SDSS magnitudes have quoted errors of less than 
±0.01 mag. Secondly, we observe that the vast major- 
ity of our global and component, broadband SEDs are 
well-matc h ed by one of the 27 spectral templates from 
[Poggiant] ([1997[ ). and thus adopt a maximum uncer- 
tainty of ±1 Hubble morphological type. This translates 
to K-correction and e-correction error contributions to 
the absolute {u—r) colors of no more than ±0.03 mag and 
±0.14 mag respectively. A further error of ±0.09 mag is 
contributed by u ncertaintes in the in ternal dust correc- 
tions derived by ([Driver et al.l[2007b[ ) and the measure- 
ment of disc inclinations. Hence, the total uncertainty 
on our absolute (u — r) colors is roughly ±0.17 mag. Of 
course, the uncertainties in our 'less confident' compo- 
nent bulge and disc colors will be significantly larger, so 
we remove these systems from our analysis of component 
colors in Section [373l 



2.6. The Impact of Unmodeled Structural Features 

The presence of unmodeled structural features in a 
number of our two-component systems contributes an 
additional source of uncertainty in the measurement of 
bulge and disc colors which is not accounted for by our 
limits on size used to identify galaxies with 'confident' 
bulge Sersic in dices and colors. Specifically, a number 
of authors ('e.g. [Laurikainen et all [20071 : [Gadottil[2008at 
iWeinzirl et al.ll2008D have recently demonstrated that ne- 
glecting to account for bars during bulge-disc decompo- 
sitions can lead to large systematic biases: bulge Sersic 
indices, scale-sizes, and fiuxes are over-estimated as the 
model bulge component is forced wider to account for bar 
flux, and the corresponding disc profile becomes steeper 
to compensate (Gadotti 20 08ah. From asample of 186 
bright nearby galaxies, lEskridge et al.l ([2000f ) estimate 
(via visual classification) that over half of all bright spi- 
rals are strongly-barred in the near-IR, although this 
fraction drops to roughly a third in the optical. Hence, 
bars are clearly an important feature to consider during 
structural decompositions. 

As a demonstration of the impact of bars on Allen 
et al.'s (2006) bulge-disc decomposition results we de- 
rive alternative fits to a number of (visually-classified) 
strongly-barred galaxies in our two-component sam- 
ple via multi-component (Sersic bulge plus Sersic bar 
lus exponential disc) fits with GALFIT version 2.0.3c 
Peng et al] [2002f ). We note that while the truncated 
Sersi c profiles offered by the rival BUDDA so ftware pack- 
age ([de Souza. Gadotti. fc dos AniosI [2004[ ) potentially 
enable a more faithful representation of bar light profiles 
(and often provided the most satisfactory fits to MGC 
galaxies in informal tests), only GALFIT currently allows 
the input of user-generated PSF functions necessary for 
fair comparisons with Allen et al.'s (2006) GIM2D analy- 
sis. The real (MGC) and model (GIM2D and GALFIT), 
B-band images and surface brightness profiles of three 
example galaxies in which the addition of a bar compo- 
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nent offered a superior fit are shown in Fig. [7] in the Ap- 
pendix, fn each case the multi-component model fit re- 
turned a much lower bulge Sersic index and B/T flux ra- 
tio than in the original, two-component model fit. Thus, 
wc recommend future, automated, bulge-disc decompo- 
sition studies incorporate the modeling of bars from the 
onset, despite the added complexity thereby introduced 
into quality control and interpretation of the output. We 
intend to explore the nature of bars in the MGC further 
in an upcoming paper (Cameron et al., in prep.). 

Finally, we note that late-type discs with bright knots 
of active star-formation also create a class of problematic 
fits. Examples of the real (MGC) and model (GIM2D 
and GALFIT), B-band images and surface brightness 
profiles of three such galaxies are displayed in Fig. [8] in 
the Appendix. A comparison is again made to alterna- 
tive fits using GALFIT, but this time with only a two- 
component (Sersic bulge plus exponential disc) model as 
these objects appear unbarred. Despite the use of iden- 
tical galaxy and PSF models, GALFIT and GIM2D re- 
turn significantly different bulge structural parameters 
in two of these examples — a refiection of their differing 
fit algorithms. The presence of complex, irregular struc- 
ture precludes clean fits with smooth elliptical isophote 
models, and identification of the global minimum be- 
comes difficult. In the examples shown, neither software 
package offers satisfactory output as each has accounted 
for some intermediate radius star-formation flux via the 
bulge model. The up-coming version of GALFIT (v3.0) 
will introduce features for mapping asymmetric, complex 
structure within the context of well-defined radial surface 
brightness profiles, and thus offers an exciting possibility 
for improved 2D bulge-disc decomposition of such galax- 
ies. 

Importantly, however, extensive visual inspection of 
galaxies in each of our key structural classes (i.e., ellip- 
tical, bulge-plus-disc, and disc-only) indicates that the 
presence of unmodeled structural components very rarely 
results in a clear mis-classification. Hence, we conclude 
that, although component bulge Sersic indices and B/T 
fiux ratios (and hence colors) for some fraction of our 
galaxies may be systematically biased, the presence of 
unmodeled bars and bright knots of active star- formation 
in the MGC i3-band images will not affect our investiga- 
tion of the role of galaxy structural type in shaping the 
global color-concentration bimodality. 

3. RESULTS 

3.1. The Global Color-Concentration Distributions of 
One- and Two- Component Galaxies 

The global color-concentration distribution of bright 
(Ms -51og/i7o < -18 mag), nearby (0.013 < z < 0.18), 
one-component galaxies in the MGC is presented in Fig. 
[21 This structural class displays strong bimodality, sepa- 
rating naturally in this parameter space into two distinct 
sub-populations — a 'red, highly-concentrated' peak and 
a 'blue, diffuse' peak — composed of elliptical and disc- 
only systems respectively. These peaks are evidently 
the foundations of those observed in the global color- 
concentrati on distribution of t he full MGC population 
(Fig. 15 in [Driver et al.l [20061 ) . although their separa- 
tion is more clearly defined in the absence of the two- 
component systems. We quantify this improvement in 
the apparent bimodality of the distribution by comparing 
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Fig. 2. — The global color-Sersic index distribution of bright 
(Mb - 51ogh70 < -18 mag), nearby (0.013 < z < 0.18), one- 
component galaxies in the MGC. Disc-only systems are represented 
by filled, blue circles and elliptical systems by filled, red circles with 
the subdivision at n = 1.5 used to define these structural types in- 
dicated by the vertical, long-dashed, black line. The corresponding 
subdivision in color at {u — r)g = 2.1 mag is indicated by the hori- 
zontal, long-dashed, black line. The volume-density distribution of 
the one-component population is revealed by the solid, black con- 
tours marked at log 4> = —5.25, —5 and —4.75 h^^ Mpc""^ per bin 
with the parameter space shown divided into a 48x48 square grid. 
The median reddening correction for internal dust extinction ap- 
plied to the disc-only systems is indicated with a black arrow. The 
mean uncertainties in the measurements of global colors and Sersic 
indices for galaxies in our sample are indicated by the magenta 
error bar at the bottom of this plot. Normalized histograms of the 
distributions of elliptical and disc-only galaxies in each variable are 
contained in the sidepanels. 

the reduced-x^ values resulting from least-squares fitting 
of pairs of bivariate Gaussian functions to the observed 
global color-concentration distributions of one- and two- 
component galaxies jv ~ 14.2) and one-component 
galaxies only {x^/u = 2.1)®. The centroids and widths 
of each peak, as derived by this bivariate Gaussian fitting 
approach, are compiled in Tabled] for reference. 

Inspection of Fig. [2] confirms that the subdivision of 
the one-component population into elliptical and disc - 
only systems at n = 1.5 adopted bv lAllen et al.l (|2006l ) 
is indeed consistent with a local minimum in this pa- 
rameter. An equivalent local minimum in color exists at 
roughly (u — r)g = 2.1 mag. Relative to these subdi- 
visions (indicated in Fig. [2]), luminous, 'bulge- less, red 
discs' and 'disc-less, blue bulges' (i.e., blue spheroids) 
are rare; specifically, we recover total volume-densities of 
only 1.7 ±0.3 and 1.1 ±0.1 x 10"'' h^g Mpc'^ respec- 
tively in the MGC. Implications of the paucity of these 
systems for galaxy formation scenarios are discussed in 
Section m 

The global color-concentration distribution of bright 

* Of course, neither of these reduced-x^ values represents a tech- 
nically good fit as expected given that both peaks are slightly as- 
symetric in shape, indicating intrinsic deviations from simple bi- 
variate Gaussian distributions. We also note that a 2D K-S test 
comparing the one and two component galaxy color-concentration 
distributions returned a probabilty of <d% that these samples were 
drawn from the same population. 
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TABLE 1 

Centroids and Standard Deviations of the One-Component 
Galaxy Color-Concentration Distribution Peaks 





(n> 




((«-r-)9> 
(mag) 


(mag) 


Elliptical 
Disc-only 


3.4 ± 0.2 
0.71 ± 0.02 


0.15 ±0.01 
0.16 ±0.01 


2.62 ±0.09 
1.44 ±0.04 


0.27 ±0.01 
0.28 ±0.04 



Note. — Values quoted are relative to the absolute magni- 
tude limit of the sample {Mg — 5 log fer o < — 18 ma^), rec alling 
that mean red and blue sequence colors (Baldrv ct al. 2004), and 
elliptical galaxy Sersic indices ([ Graham fc Guzman 2003), are lu- 
minosity dependent. 



(Mb -51og/i7o < -18 mag), nearby (0.013 <z < 0.18), 
two-component galaxies is presented in Fig. [3) These 
bulge-plus-disc systems constitute an intermediate color- 
concentration class, spanning the two peaks of the one- 
component galaxy bi modality. As mentioned earlier, 
iDrorv fc Fished ()2007f ) identified a relationship between 
visually-classified bulge type (i.e., 'pseudobulge' or 'clas- 
sical' bulge)^ and position in the color-concentration 
plane. Their sample of 39 SO-Sbc galaxies is overlaid 
against the MGC distribution in Fig. [3] for comparison, 
and a close agreement is evident. (Note that we ad- 
just the colors quoted in Drory & Fisher's (2007) Ta- 
ble 1 for the effects of internal dust reddening, and for 
an offset of ~0.2 mag between their integrated Sersic 
profile colors a n d the SDSS model colors used here.) 
iDrorv fc Fished ()2007D reveal that Sabc, pseudobulge 
galaxies may be isolated from SO/Sabc, classical bulge 
(and SO pseudobulge) galaxies in this parameter space 
using cuts dX n — 1.5 and {u — r)g — 2.2 mag (~2.1 
mag after the relevant corrections). Motivated by these 
results we subdivide our two-component population into 
two structural types, 'bridging' and 'red peak', defined 
via a cut at {u — r)g — 3.22 — 2.75 log n (indicated in 
Fig. [3]). Visual inspection of the i?-band MGC images 
confirms these structural types correspond to morpho- 
logically l ate-type and early-type disc respectively, as re- 
ported bv IDrorv fc Fished (|2007f) . 

Examining the two-component, red peak systems we 
note that their global colors are significantly bluer than 
those of the one-component, red peak (i.e., elliptical) sys- 
tems, which appears to indicate a younger mean stellar 
population age for these galaxies (modulo the effects of 
metallicity) . However, this offset is almost entirely in- 
troduced by our correction for dust reddening, which (as 
noted earlier) may not be appropriate for these early- 
type discs. Hence, we caution against placing excessive 
emphasis on this color offset in the interpretation of our 
results. We also note that the two-component, red peak 
systems display highly-concentrated light distributions 
similar to those of the ellipticals. Although this appears 
to argue for the role of violent-relaxation processes dur- 
ing their formation, the known prevalence of (secular 
evolution -built) pseudobulg es amongst earl y-type disc 
galaxies (lErwin et al.l l2003t iLaurikainen. S alo. fc Butal 

^ Although visual b ulge morpholog y has been shown to correlate 
with V — H c olor ( Car oUo et al.l[2002l') and mid-infrared bulge color 
lIFisheii I2OO6 I wc would caution that dynamical data (indicating 
the degree to which these structures are rotationally-supported or 
otherwise) would allow for more confident classification of bulge 
type. 
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Fig. 3. — The global color-Sersic index distribution of bright 
[Mb - 51og/i7o < -18 mag), nearby (0.013 < z < 0.18), two- 
component galaxies in the MGC. These bulge-plus-disc systems 
are represented by filled, orange circles. For reference, the volume- 
density distribution of the one-component population is revealed by 
the solid, black contours marked at log <j> = —5.25, —5 and —4.75 
h^Q Mpc~^ per bin with the parameter space shown divided into a 
48x48 square grid. The median reddening correction for internal 
dust extinction applied to these bulge-plus-disc systems is indicated 
with a black arrow. Drory & Fisher's (2007) sample of 39 SO-Sc 
galaxies (with colors adjusted for dus t red dening and a magnitude 
system offset, as described in Section |3. Il l is overlaid for compari- 
son. Their (visually-classified) pseudobulge galaxies are marked as 
green squares and their classical bulge galaxies as purple triangles 
with open symbols denoting SOs and solid symbols Sa-Sc types. 
Our subdivision of the two-component population into 'bridging' 
and 'red peak' structural types at (it — r)g = 3.22 — 2.751og(n) is 
indicated by the long-dashed, black line. For comparion, Drory & 
Fisher's (2007) subdivisions at n = 2.5 and (u — r)g = 2.2 mag (ad- 
justed to 2.1 mag) are indicated by the short-dashed, black lines. 
The mean uncertainties in the measurements of global colors and 
Sersic indices for galaxies in our sample are indicated by the ma- 
genta error bar at the bottom of this plot. Normalized histograms 
of the two-component galaxy distribution in each variable are con- 
tained in the sidepanels. 



120051: IDrorv fc Fished l2007f ) challenges this interpreta- 
tion. 

The two-component, bridging population systems are 
generally similar in color to their one-component, disc- 
only couterparts, albeit slightly redder, whilst their 
global Sersic indices are significantly higher. These re- 
sults are broadly consistent with the interpretation of 
these galaxies as pseudobulge systems given that N-body 
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Fig. 4. — The B-band luminosity functions of key structural 
types identified in the MGC. Namely, one-component {light green) 
and two-component (orange) systems, and their corresponding 
sub-types: elliptical (red) and disc-only (blu e); bridging {dark 
green) and red peak {purple) (see Section 13.11 for definition). In 
each case, the measured number density of galaxies in each 0.5 
mag bin is indicated by a solid square with Poisson error bars, 
and the corresponding best-fit Schechter luminosity function by a 
solid line. The open squares with error bars and dotted lines illus- 
trate the effect of removing the internal dust extinction correction 
applied to the two-component, red peak systems. The Schechter 
function parameters of each fit are contained in Table [2l 

simulations demonstrate the effectiveness of secular evo- 
lution in building up comparable central mass conccn- 
trations over 1-2 Gyr t imescales (e.g.^thanassQula 2005.; 
iDebattista et al.ll2006f ). We will return to the interpre- 
tation of our results in light of galaxy formation theory 
in Section 0] — but first we recover luminosity functions 
and total stellar mass densities for each structural type 
(Section 13. 2p . and explore the connection between bulge 
and disc colors in two-component systems (Section I3.3|) . 

3.2. Luminosity Functions and Stellar Mass Densities 
of Key Structural Types 

The luminosity functions of the key MGC structural 
types identified in Section 13. H — namely, one-component 
(elliptical and disc-only) and two-component (bridging 
and red peak) — are presented in Fig.|4l Each of these was 
modeled using a Schechter luminosity function with the 
best-fit Schechter parameters and integrated luminosity 
densities contained in Tabled) 

Inspection of the luminosity functions by structural 
type in Fig. [4] reveals that two-component galaxies domi- 
nate (in number density) over one-component galaxies at 
the brightest luminosities {Mb — 51og/i7o < —20 mag), 
whereas one-component galaxies dominate at intermedi- 
ate luminosities (—20 < Mb — 51og/i7o < —18 mag). 
Furthermore, over the luminosity range explored, disc- 
only systems are more common than ellipticals, and red 
peak systems are more common than bridging systems.^'' 
No previous study has examined luminosity functions by 

However, if the dust attenuation corection is removed from the 
red peak systems (as discussed in Section [2.4|l . bridging systems are 
more common at the brightest luminosities. 



structural type using an equivalent bulge-disc decompo- 
sition approach. However, we note that our results are 
qualitatively consistent with previous luminosity func- 
tions by visual morphological class (given the expected 
morphological composition of each structural type). In 
particular, the B-band luminosity functions compiled by 
Ide Lapparentj (^003) reveal that Sa/Sb (two-component, 
red peak and bridging) galaxies dominate at the bright- 
est luminosities, while Sc (one-component, disc-only; and 
two-component, bridging) and Sd/Irr (one-component, 
disc-only) galaxies dominate at intermediate and faint 
luminosities respectively. 

We also note that the total luminosity density of all 
structural types derived here, Jb — 2.6 ± 0.3 x 10^ hjo 
Lq Mpc~^ is similar to that for all structural component s 
(i.e., bulges and discs) derived bv [Driver et al.l ()2007b[ ). 
Jb = 3.0 ± 0.2 X 10* hjo Lq Mpc'^. Both studies 
utilize the MGC structural catalogue, however, the to- 
tal luminosity density obtained here is slightly smaller 
due largely to the reclassification of B/T > 0.6 galax- 
ies as ellipticals (intended to remove 'false discs'; see 
Section I2.1.ip — when returned to the elliptical popula- 
tion these galaxies are no longer dust corrected, and 
hence become less luminous. We therefore also note 
the revised luminosity density of elliptical galaxies de- 
rived here, Jb = 0.27 ± 0.07 x 10* hrg Lp, Mpc~^ whic h 
is higher than that recovered bv [Driver et al.l ()2007b[ ). 
Jb = 0.14 ± 0.02 x 10* hjo Lq Mpc'^. 

We also compute stellar mass densities using the model 
relationships betwe en (m — r) color and mass-to-light ra- 
tios given by Bell et al.l (j2003f ). We thereby recover total 
stellar mass densities in each structural type of 0.7 ± 0.1, 
1.3±0.1, 1.8±0.2 and l.liO.l xlO* /170 Mq Mpc'^ for el- 
liptical, disc-only, red peak and bridging systems respec- 
tively. Again we note that our total stellar mass density 
for aU structural types, 4.9 ± 0.3 x lO* hjp Mq Mp c"^, 
is similar to that recovered by iDriver et al.l ()2007bf ) for 
all structural components (i.e., bulges and discs) in the 
MGC, 5.2 ±0.4 xlO* /i7o Mq Mpc'^. 

The comparable luminosity and mass densities of one- 
and two-component systems revealed here highlights the 
importance of considering the fundamental distinction 
between elliptical, disc-only and bulge-plus-disc systems 
in order to properly understand galaxy formation and 
evolutionary processes. These luminosity functions by 
structural type, and those by sep a rate str uctural com- 
ponent presented in iDriver et al.l (|2007bl ). offer valu- 
able constraints for galaxy formation models, and should 
serve as key reference points for semi-analytical simula- 
tions. 

3.3. Relationship between Bulge and Disc Colors 

In Section [01 we demonstrated that the dichotomous 
disc-only and elliptical (or 'bulge-only') galaxy popula- 
tions are the foundations of the bulge, low-n and red, 
high-n peaks of the global color-concentration bimodal- 
ity, and that bulge-plus-disc galaxies span these two 
peaks — as expected given the correlation between global 
color-conce ntration and e yeball morphological type re- 
por ted bv IDriver et al.l j2006D . However, as noted 
by iDrorv fc Fished (|2007D . Driver et al.'s (2006) hy- 
pothesis that the global properties of bulge-plus-disc 
systems are the result of mixing red, high- 71 bulges 
and blue, low-n discs in varying degrees is contra- 
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TABLE 2 

Luminosity Function Parameters for Key Structural Types in the MGC 





Mg - 5 log hjo 


a 




is 


P 




(mag) 




(10-3 h^Q Mpc-3 [0.5 mag]-i) 


(lOS hro Lq Mpc-3) 


(10« hro Mq Mpc-3) 


Ono-Component (All) 


-19.05 ± 0.07 


-1.03 ±0.04 


5.8 ±0.2 


1.3 ± 0.1 


2.0 ±0.2 


Disc-only 


-19.14 ±0.07 


-1.08 ±0.04 


4.0 ±0.2 


1.0 ±0.1 


1.3 ±0.1 


Elliptical 


-18.93 ± 0.08 


-0.69 ±0.04 


1.5±0.1 


0.27 ± 0.07 


0.7 ±0.1 


Two-Component (All) 


-19.82 ± 0.08 


-0.61 ±0.04 


3.3 ±0.2 


1.3 ±0.1 


2.9 ±0.3 


Red peak 


-19.69 ± 0.09 


-0.20 ±0.05 


1.6 ±0.1 


0.59 ± 0.06 


1.8 ±0.2 


Red peak (non-DC) 


-19.10 ±0.08 


-0.39 ±0.05 


1.8 ±0.1 


0.37 ±0.04 


1.4 ±0.2 


Bridging 


-19.59 ± 0.08 


-0.84 ±0.04 


2.0 ±0.2 


0.68 ± 0.07 


1.1 ±0.1 



Note. — All quoted uncertainties are the random errors in the Schechter function fits. Based on mock 2dFGRS NG P catalogues 
llCole et al.lll99^ . we estimate that the potential systematic errors on these values due to cosmic variance amounts to 13% ijDriver et al.l 
120051) . A similar level of uncertainty exists in the elliptical, two-component (all) and red peak param eters d ue to the potential for mis- 
classification of ellipticals as bulge-plus-disc systems (see the discussion regarding 'false discs' in Section 12.1.11 and Fig. |6ll. 



dieted by t he close relationship between bulge and 
disc colors (jPeletier fc Balcellsl Il996t iMacArthur et all 
[2004). Here we confirm this relationship for the bulge 
and disc colors of MGC two-component galaxies (with 
'confident' colors), as shown in Fig. [5l In Peletier 
& Balcells' (1996) study of 30 highly-inchned, SO- 
Sbc galaxies the authors measured component colors 
via wedge-shaped apertures designed to avoid biases 
due to dust attentuation in the plane of each disc. 
They thereby measured a ('dust-free') mean bulge-disc 
color offset of (([/ - R)b - (U - R)d) = 0.126 mag, or 
{{B — R)i, — [B — R)d) — 0.045 mag in a n alternative fil - 
ter combination. More recently, MacArth ur et al.| (j2004l ) 
reported a substantially larger mean bulge-disc color dif- 
ference of {{B — R)b — {B — R)d) = 0.29 mag from an 
investigation of color gradients in a sample of 172 low- 
inclination, SO-Sd/Irr galaxies. In that study component 
colors were estimated from the averag e colors at radii of 
0-0. 5/i for bulges and 1.5-2.5/i for discs. IMacArthur et al] 
(|2004 ) speculated that their disagreement with Peletier 
& Balcells' (1996) offset was either due to internal dust 
extinction, or a bias in the wedge-shaped aperture colors 
of the latter due to vertical gradients in stellar popula- 
tion age and metallicity within discs. Here we recover an 
offset of {{u — r)b — {u — r)d) = 0.22 ± 0.02 mag using 
our dust attenuation corrected ([Driver et al.ll2008h col- 
ors, which is indeed also larger than Peletier & Balcells' 
(1996) offset^^. Removing the dust correction from our 
red peak systems results in an even larger offset for that 
population of {{u — r)b — {u — r)d) = 0.27 ± 0.04 mag. 

4. CONCLUSIONS AND DISCUSSION 

Our key conclusions may be summarized as follows. 
Via bulge-disc decomposition and quantitati ve morpho- 
logical classification (i.e., the logical filter in I Allen et al.l 
[2006) of luminous {Mb — bloghm < —18 mag), nearby 
(0.013 < z < 0.18) galaxies in the MGC, we have demon- 
strated that: 

(i) one-component systems (predominantly elliptical 
and disc-only galaxies) display strong bimodality in the 
color-concentration plane, contributing to a red, high 
Sersic index peak and a blue, low Sersic index peak re- 
spectively (with total stellar mass densities of 0.7 ± 0.1 
and 1.3 ± 0.1 x 10* hm Mq Mpc'^ in each); 

We suppose a close agreement between (u — r) and {U — R) 
color differences based on the relevant filter transformation equa- 
tions presented in Blanton & Roweis' (2007) Table 2. 
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Fig. 5. — The relationship between the bulge and disc colors 
of two-component galaxies in the MGC. Only systems with 'con- 
fident' component colors are shown on this plot. Each bulge is 
marked by a circle, colored dark green if the bulge is embedded 
in a bridging population galaxy and purple if the bulge is embed- 
ded in a red peak population galaxy. The estimated uncertainties 
in the measurements of bulge and disc colors for galaxies in our 
sample are indicated by the magenta error bar. 

(ii) luminous, 'bulge-less, red discs' and 'disc-less, blue 
bulges' (i.e., blue ellipticals) are rare; (with volume- 
densities of 1.7 ± 0.3 and 1.1 ± 0.1 x lO"'^ h^g Mpc'^ 
respectively); 

(ill) two-component systems (bulge-plus-disc galax- 
ies) constitute a population of intermediate color- 
concentration, spanning the two peaks of the bimodality 
defined by one-component systems; 

(iv) component bulge and disc colors are tightly cor- 
related (with mean offset {{u — t)i, — {u — r)d) = 0.22 ± 
0.02 mag); and 

(v) the two-component systems may be subdivided in 
a physically-motivated manner at {u — r)g < 3.22 — 
2.75 log n into bridging (late- type disc) and red peak 
(early-type disc) populations (with total stellar mass 
densities of 1.1 ± 0.1 and 1.8 ± 0.2 x 10* hro Mq Mpc'^ 
in each). 

We now investigate how these results compare against 
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other recent studies of galaxy structure, and discuss their 
interpretation within the hierarchical clustering scenario 
of galaxy formation, focussing on elliptical, disc-only, 
and bulge-plus-disc systems in turn. Finally, we out- 
line future work, and discuss the challenges facing future 
low resolution, automated structural decomposition cam- 
paigns. 

Ellipticals — The red colors and highly centrally- 
concentrated light profiles of luminous ellipticals in the 
MGC are consistent with many other observational stud- 
ies, including a large sample of elliptical galaxies in the 
SDSS selected via a quantitative m orphological approach 
and examined by iBernardi et al.l (j^OOSa.b c d). More- 
over, these properties are broadly consistent with the 
modern picture of hierarchical formation in which bright 
ellipticals are assembled at late e pochs via 'dry merg- 
ers' with other early-type sy stems ()Khochfar fc BurkertI 
I2003t iDe Lucia etall |2006j. Only their progenitors 
(i.e., lower luminosity/mass ellipticals) are built up early 
via similar-size d disc mergers as in the original hier - 
archical model (jToomre fc ToomrdfTOTl lBarne£lll992l ). 
The modern scenario is supported by direct observa- 
tions of early-type merg ers at intermediate redshifts 



tions 01 early-type mergers at mte rmediate redsnitts 
(|van Dokkum et al.|[T999t " lTran et al.ll2005f ). and N-body 



simulations confirming that 'dry mergers' can repro- 
duce the st ructural an d kinematic properties of nearby 
ellipticals (|Naab fc Bu rkcrt 2 0031) and preserve the 
Fundamental Plane (iNipoti. LMulrillo. fc Ciotlil l2003l : 
iBovlan-Kolchin. Ma. fc Quataert]l2005( l. 

The 'dry merger' scenario is also purported to account 
for the observ ed build up of mass on the red sequence 
since z ~ 1 ijBell et al.l [200l iBlantonI [2001 . How- 
ever, studies showing minimal evolut ion of the early -type 
stellar mass-size relation (.Mcintosh et al.l l2005l) — not 



prese rved by 'dry merging' (|Nipoti. Londrillo. fc Ciottil 
|2003| ) — and constraints on the frequency of recent merg- 
ing in l uminous galaxies from p air counts and asymmetry 
indices (|de Propris et al.|l2007^ challenge this hypothesis, 
suggesting a greater role for (non-merger) environmental 
processes in the recent evolution of the red sequence. 

The study of galaxy structure via automated bulge- 
disc decomposition of large, intermediate-to- high redshift 
datasets, followed by comparison against local samples 
such as the MGC, may enable the degeneracy between 
early-type, bulge-plus-disc (SO) and predominantly el- 
liptical (E) systems implicit in studies of the 'red se- 
quence' to be broken — thereby providing more accurate 
constraints on elliptical galaxy evolution. However, the 
handling of 'false disc' candidates (see Section r2.1.ip re- 
quires improvement. In the meantime, the MGC ellip- 
tical population should serve as a key local benchmark, 
which semi-analytical models must be able to reproduce 
whilst simultaneously accounting for the properties of 
disc-only and bulge-plus-disc galaxies as well. 

Disc-only systems — The blue colors and diffuse (i.e., 
low-n) light profiles of predominantly disc-only systems 
in the MGC are also consistent with previous obser- 
vations of disc-only galaxies, including a large sample 
of bulge-less, e dge-o n discs identified in the SDSS by 
iKautsch et al.l ()2006f l. These properties are also quali- 
tatively consistent with expectations from the hierarchi- 
cal formation scenario in which disc-dominated galax- 
ies form early, but experience a passive evolutionary his- 



tory (iFall fc Efstathioul[l980l : IMo. Mao, fc Whitelll998fl . 

However, the prevalence of disc-only systems in the lo- 
cal universe is difficult to reconcile with the expected 
merger history of galaxies in CDM cosmologies. Specifi- 
cally, in our analysis of the MGC data we identify disc- 
only systems as contributing '^30% of the total stel- 
lar mass in nearby disc galaxi es, which is compara - 
ble to the number reported by IKautsch et al.l (|2006D . 
In a theoretical study of halo merger historie s in a 
ACDM universe, iKoda. Milosavlievic. fc Shapirol ()2007f ) 
reveal that the observed frequency of bulgeless discs can 
only be reconciled within the hierarhical formation con- 
text provided spheroid production only occurs in merg- 
ers involving mass ratios greater than 0.3. Moreover, 
'Stew art et al.l ()2008f ) reveal equivalent limits on the fre- 



quency of minor mergers, which are expected to fuel cen- 
tral st ar-formation, leading to the creation of pseudob- 
ulges (|Herauist fc Mihoslll995l ). 

One may question whether these disc-only galaxies 
identified in the MGC are indeed entirely 'bulge-less', 
or whether they simply possess very faint 'bulges', in- 
significantly brighter in central surface brightness than 
the discs in which they are embedded. This is indeed 
possible given the li mitations of our photometric de- 
composition pipeline ([Allen et al.ll2006l ). and we further 
note the po tential for dust extinction to obscur e com- 
pact bulges ([Tuffs et al.ll2004l : [Driver et al.ll2007bl ). Nev- 
ertheless, even if a population of very faint bulges is hid- 
ing within our disc-only population, their properties are 
vastly different to those of the bright, high surface bright- 
ness classic al bulges predicted by t he hierarhical cluster- 
ing model ([Baugh. Cole, fc FrenkI [l996l ) . We also note 
that i n a high resolution study of 19 nearby late-type 
discs, iBoker. Stanek. fc van der Marell ([2003) observed 
that most of their sample could be well-fit using a sin- 
gle Sersic profile model, and the vast majority of central 
flux excesses detected could not be attributed to 'what 
one woul d generally consider to be a bul ge component'. 
Similarlv. lWvse. Gilmore. fc Franxl ([19971 ) point out that 
of the four largest disc galaxies in the local group, two 
(M33 and the LMC) do not possess convincing classical 
bulges. 

However, it is certainly important to account for obser- 
vational limitations when comparing output from semi- 
analytical simulations against empirical datasets in or- 
der to derive robust conclusions. One example of how 
thi s may be achieved is the analysis pipeline developed 
bv iHaussler et al.l ([2007f ) to generate reahstic simulated 
galaxy images from model predictions, which may be 
used to compute 'observed' structural parameters. 

Bulge-plus- disc systems — The close relationship be- 
tween the bulge and disc colors of two-co mponent 
galaxi e s — as well as recen t low resolution ([Gadottil 
I2008U 'Weinzirl et al.' '20081 and hi gh resolution s tud- 
ies ([Drory fc Fisher 2007; Laurikai nen et all [20071 ) re- 
vealing the disc-like nature of many blue and red 
bulges — points towards a scenario in which both sec- 
ular evolution and merging play key roles in bulge 
formation, in contrast to the c o nvent ional hierarchi- 
cal view ([Baugh. Cole, fc FrenkI I1996D. Recerit N- 
body s imulations (Bournaud. Combes, fc Semelinll2005l : 
[Athanassoula 2005; Debattista et al.„ .200(31 ) certainly in- 
dicate that secular evolutionary processes are efficient 
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at establishing pseudobulge-like, central mass concen- 
trations. However, the distinction between classical 
bulge and pseudobulge systems may not be as clear 
as is often suggested — bas ed upon SAURON o bserva- 
tions of 24 Sa/Sab galaxies. iPeletier et alJ (|2008^ report 
the co-existence of both an old, velocity-dispersion sup- 
ported 'classical' component and a young, rotationally- 
supported 'pseudobulge' component within the nuclei of 
these early-type discs. 

High redshift observations offer some further impor- 
tant clues as to the formation pathways of bulge-plus- 
disc galaxies, although cosmic variance contributes a 
large source of uncertainty to the results. For instance, 
from the struct ural decomposition of 24 8 galax ies at 
0.1 < z < 1.3, iDomfnguez-Palmero fc Balcellsl (|2008[ ) 
reveal that the correlation betweeen bulge and global 
color is in place by intermediate redshifts. They also 
identify 60% of their bulges as displaying red colors con- 
sistent with passively-evolving stellar populations, and 
40% as displaying bluer colors indicative of more recent 
star-formation. They speculate that the formation of 
red bulges coincided with a truncation of galaxy star- 
formation that did not destroy th e disc. Indeed, the role 
of black hole and AGN formation (jSomerville et al.ll2008f) 
cannot be overlooked in any complete model of bulge for- 
mation given the relationshi ps bet w een bulge and black 
hole properties jM agorrianj Il998t iFerrarese fc MerrittI 
2000HGebhardt et al...200Q,). In contrast to the results o f 
Dommguez-Palmero fc Balcellsl (|2008D . iKoo et al] (|2005l ) 
report a much higher fraction (85%) of red bulges at high 
redshift in their sample of 71 cluster and 21 field galaxies 
at z ~ 0.8, more consistent with a merger-built origin. 
However, the observed fraction is perhaps influenced by 
their selection of luminous bulges (rather than simply 
globally luminous galaxies). 

Future work, and challenges for low resolution, au- 
tomated structural decomposition studies — Semi- 
anal ytical models of hierarchical galaxy formation 
(e.g. Ivan Kampen. Jimenez, fc Peacocklll999l : ICole et al.l 
now have the power to trace the formation and 
evolution of both bulge and disc components separately 
(| Almeida. Baugh. fc Lacevll2007[ ). Moreover, wide-field 
galaxy surveys (e.g. SDSS, MGC, GEMS) now achieve 
sufficient survey volumes and imaging resolutions 
to facilitate large-scale, autom ated/semi-automated , 
bulge-disc deconiposition s tudies fTasca fc White"2005'; 
lAllen et all [200l iGadottil [2008b: Weinzirl et al. 2008). 



A number of recent studies also offer important — albiet 
relatively small sample — benchmarks of galax y structure 
in th e intermediate-to-high redshift universe ()Koo et al.l 
I2005t IDomfnguez-Palmero fc Balcellsl I2008D . Hence, 
there currently exist unprecedented opportunities to 
unify theory and observation via detailed statistical 
studies of galaxies by structural type and component. 

However, many observational challenges in this en- 
deavor still remain. Advanced, automated, data pro- 
cessing and quality control procedures are required to 
account for the impact of secondary structural features 
(e.g. bars and/or bright knots of active star-formation) 
during bulge-disc decompositions whilst ensuring mean- 
ingful output. Automation of these processes is manda- 
tory due to the vast sample sizes necessary to successfully 
limit cosmic variance uncertainties, provide accurate sta- 
tistical constraints on the relevant structural parameters, 
and enable investigation of the role of environment and 
mass. Finally, we note that a meaningful comparison 
between theory and observation will also require care- 
ful handling of the relevant observational selection bi- 
ases and measurement errors, and the use of quantitative 
classification techniques for defining equivalent morpho- 
logical sub-types. Fortunately, quality, wide-field, near- 
IR survey data is becoming increasingly available (e.g. 
UKIDSS, GAMA), which should help to hmit the im- 
pact of irregular star-formation features and dust atten- 
uation. Furthermore, the next generation of bulge-disc 
decomposition software is expected to introduce many 
new capabilities (as anticipated in the up-coming GAL- 
FIT v3.0). 



The Millennium Galaxy Catalogue consists of imaging 
data from the Isaac Newton Telescope and spectroscopic 
data from the Anglo Australian Telescope, the ANU 2.3- 
m, the ESO New Technology Telescope, the Telesco- 
pic Nazionale Galileo and the Gemini North Telescope. 
The survey has been supported through grants from 
the Particle Physics and Astronomy Research Coun- 
cil (UK) and the Australian Research Council (AUS). 
The data and data products are publicly available from 
http://www.eso.org/jliske/mgc/ or on request from J.L. 
or S.P.D. E.G. acknowledges partial financial support 
from the Australian Research Council Discovery Project 
Grant DP0451426. 
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APPENDIX 

This Appendix contains real (MGC) and model (GIM2D/GALFIT), i?-band images and surface brightness profiles 
of example galaxies illustrating three problematic fit scenarios — false discs (Fig. |6]), strong bars (Fig. [7|), and bright 
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knots of active star-formation (Fig. [8]). 
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Fig. 6. — Comparison of the GIM2D, two-component (Sersic bulge plus exponential disc) and one-component (Scrsic-only) model fits to 
three bright MGC galaxies identified as two-component systems by the logical filter, but believed to be ellipticals based on visual inspection 
of their B-band images and surface brightness profiles. In each case the elliptical isophote, semi-major axis, surface brightness profile of 
the real galaxy (recovered using the Starlink ELLPRO package) is indicated by open, blue squares and 3cr error bars. The corresponding 
2D model surface brightness profile (recovered from the PSF-convolved, GIM2D model output image, again using the ELLPRO package) is 
indicated by filled, red triangles, while the broken, black lines reveal the (non-PSF-convolved) contributions of each component separately. 
Postage stamps of the B-band, MGC images, GIM2D model images and the corresponding residuals are also presented in greyscale with 
a logarithmic weighting function. In each of these examples the addition of the outer exponential component results in only a trivial 
improvement over the one-component fit, as reflected in the relative reduced-x^ values. In the first example the improvement arises from 
fitting a mild asymmetry in the light distribution at intermediate radii, perhaps created in a recent merger event. In the second and third 
examples the improvement arises from fitting a slight fiux excess at intermediate-to-large radii, perhaps due to the presence of unresolved 
satellite(s). As the two-component model is intended for 'genuine' (i.e., rotationally-supported, fiattened) exponential discs in this study, 
the Sersic-only fit was preferred for these suspected ellipticals, which were, therefore, returned to the one-component population. 
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Fig. 7. — Comparison of the GIM2D, two-component (Sersic bulge plus exponential disc) and GALFIT, three component (Sersic bulge 
plus Sersic bar plus exponential disc) model fits to three bright MGC galaxies observed to contain strong-to-intermediate bars upon visual 
inspection of their B-band images and surface brightness profiles. In these examples the two-component model fails to provide a satisfactory 
representation of these intrinsically three (or more) component systems. The presence of the bar results in systematic biases in the structural 
parameters recovered in the two-component fit. Namely, a strong increase in the bulge Sersic index, half light radius, and flux, and a weak 
decrease in the disc scale-length and disc fiux. The labour-intensive option of refitting the MGC-BRIGHT sample with multi-component 
models using GALFIT or BUDDA will be explored in a future paper (Kelvin et al., in prep.). 
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Fig. 8. — Comparison of the GIM2D and GALFIT, two-component (Sersic bulge plus exponential disc) model fits to three bright MGC 
galaxies observed to be undergoing strong, localised, clumpy star-formation within the inner-to-intermediate disc region. These examples 
illustrate the difficulty of recovering sensible bulge structural parameters from late-type disc galaxies, which typically possess only weak 
bulges and display much irregular structure within the disc. Consequently, the most physically-meaningful, two-component fit may not 
always lie at a global, reduced-x'^ minimum, or it may be difficult to identify due to many nearby reduced-x^ minima within the model 
parameter space. Typically, this results in a systematic over-estimation of bulge Sersic index, half light radius, and flux. In the first 
and third example GALFIT appears to have converged upon a better flt, and in the second example the GIM2D flt may be preferred. 
Accounting for the effects of irregular patterns of star-formation in the disc remains a challenge for contemporary, structural decomposition 
software, expected to be addressed by the up-coming GALFIT version 3.0. 



